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Abstract 

The Spherical Neutral Detector (SND) operates at VEPP-2M collider in Novosibirsk 
studying e + e~ annihilation in the energy range up to 1.4 GeV. Detector consists of 
a fine granulated spherical scintillation calorimeter with 1632 Nal(Tl) crystals, two 
cylindrical drift chambers with 10 layers of sense wires, and a muon system made of 
streamer tubes and plastic scintillation counters. The detector design, performance, 
data acquisition and processing are described. 



1 Introduction 



For more than 25 years the VEPP-2M e + e~ collider has been operating in 
BINP, Novosibirsk, in the center-of-mass energy range 2E = 0.36 -r- 1.4 GeV 
[1]. Until now its maximum luminosity of L — 3 • 10 30 cm~ 2 s~ 1 at E = 510 
MeV was a record for the machines of this class. During this period several 
generations of detectors performed experiments at VEPP-2M. Much of the 
current data on particles properties [2] at low energy region were obtained in 
these experiments. 

1 E-mail: achasov@inp.nsk.su, FAX: +7(383-2)34-21-63 
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The SND detector is an advanced version of its predecessor - the Neutral 
Detector (ND) [3,4], which completed its five-year experimental program in 
1987 [5]. The SND detector [6] operates at VEPP-2M since 1995 [7-10]. Its 
main part is a three-layer spherical calorimeter consisting of 1632 crystals 
Nal(Tl). The SND is a general purpose detector optimized for multi-photon 
final states. 



2 General overview 

The detector layout is shown in Figs. 1, 2. Electron and positron beams collide 
inside the beryllium beam pipe with a diameter of 2 cm and 1 mm wall. The 
beam pipe is surrounded by tracking system consisting of two drift chambers 
and a cylindrical scintillation counter between them. The solid angle coverage 
of the tracking system is about 98% of 4ir. 

The three-layer spherical electromagnetic calorimeter based on Nal(Tl) crys- 
tals surrounds the tracking system. The total calorimeter thickness for par- 
ticles originating from the interaction region is 34.7 cm (13.4 Xq) of Nal(Tl) 
and the total solid angle is 90% of 47T. 

Outside the calorimeter a 12 cm thick iron absorber is placed in order to atten- 
uate the residuals of electromagnetic showers. It is surrounded by segmented 
muon system which provides both muon identification and cosmic background 
suppression . Each segment consists of two layers of streamer tubes and a plas- 
tic scintillation counter, separated from the tubes by 1 cm iron plate. The iron 
layer between the tubes and the counter reduces the probability of their si- 
multaneous firing by photons produced in e + e~ collisions to less than 1% for 
700 MeV photons. 



3 Calorimeter 

3. 1 Calorimeter layout 

Spherical shape of the SND calorimeter provides relative uniformity of re- 
sponse over the whole solid angle. Pairs of counters of the two inner layers 
with thickness of 2.9 and 4.8 X (X = 2.6 cm) are sealed in thin (0.1 mm) 
aluminum containers, fixed to an aluminum supporting hemisphere (Fig. 3). 
Behind it, the third layer of Nal(Tl) crystals, 5.7 X thick, is placed. The gap 
between the adjacent crystals of one layer is about 0.5 mm. The total number 
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of counters is 1632, the number of crystals per layer varies from 520 to 560. 
The total mass of Nal(Tl) is 3.5 t. 

The polar angle coverage of the calorimeter is 18° < 9 < 162°. The calorimeter 
is divided into two parts: "large" angles 36° < 9 < 144° and "small" angles 
- the rest. The angular dimensions of crystals are A0 = A9 = 9° at "large" 
angles and A0 = 18°, A9 = 9° at "small" angles. Each layer of the calorimeter 
consists of crystals of eight different shapes. 

The crystal widths approximately match the transverse size of an electromag- 
netic shower in Nal(Tl). Two showers can be distinguished if the angle between 
them is larger than 9°. If this angle exceeds 18° the energies of the showers 
can be measured with the same accuracy as for isolated shower. A high granu- 
larity of the calorimeter is especially useful for the detection of multi-particle 
events. For example, the detection efficiency for 7-photon events of the process 
— > r^y, rj — > 37T° is close to 15%. 

The light collection efficiency varies from 7% to 15% for crystals of different 
calorimeter layers. The scintillation light signals from the crystals are detected 
by vacuum phototriodes [11] with a photocathode diameter of 17 mm in the 
first two layers and 42 mm in the third layer. The average photocathode 
quantum efficiency is about 15% and the mean tube gain is about 10. 

The electronics of the calorimeter (Fig. 4) consists of 

(1) the charge sensitive preamplifiers (CSA) with a conversion coefficient 
of 0.7 V/pC, 

(2) 12-channel shaping amplifiers (SHA) with a remote controlled gain that 
can be set to any value in the range from to a maximum with a reso- 
lution of 1/255. 

(3) 24-channel 12-bit analog to digital converter (ADC) with a maximum 
input signal U max = 2V, 

Each calorimeter channel can be tested using a precision computer-controlled 
calibration generator. The amplitude of its signal can be set to any value 
from to 1 V with a resolution of 1/4096. The equivalent electronics noise of 
individual calorimeter channel lies within the range of 150 -r- 350 keV. 

To produce signals for the first-level trigger the SHAs are collected into 160 
modules corresponding to 160 towers (4 crystals from three layers located 
one under another). The towers divide calorimeter into 20 sectors of 18° in 
azimuthal direction and into 8 rings of 18° in polar direction. 
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3.2 Energy calibration 



Calorimeter is calibrated using cosmic muons [12] and e + e~ — > e + e~ process 
events [13]. A fast preliminary calibration based on cosmic muons gives the 
constants for calculation of energy deposition in the calorimeter crystals. These 
constants are used for levelling responses of all crystals in order to obtain 
an uniform first-level trigger energy threshold over the whole calorimeter and 
represent seed values for the more precise calibration procedure, using e + e~ — > 
e + e~ events. The cosmic calibration procedure is based on the comparison of 
the experimental and simulated energy depositions in the calorimeter crystals 
for cosmic muons. The detailed description of the method is given in [12]. 
The statistical accuracy of 1% in the calibration coefficients was achieved. 
After cosmic calibration the peak positions in the measured energy spectra 
for photons and electrons agree at a level of about 1% with the actual particle 
energies (Fig. 5). 

The cosmic calibration procedure performed weekly, between the experimen- 
tal runs takes less than 5 hours. For a one-week period between consecutive 
calibrations the coefficients stability was better than 1.5%. Changes in the 
conversion constants were also monitored daily, using the calibration genera- 
tor. They could either drift slowly with electronics gain changes or show large 
leaps due to replacement of broken electronics modules. 

To achieve highest possible energy resolution, the precise OFF-LINE calibra- 
tion procedure, based on analysis of e + e~ — > e + e~ events, was implemented. 
The calibration coefficients are obtained by the minimization of the r.m.s. of 
the total energy deposition spectrum for the electrons with fixed energy. The 
detailed description of the procedure is given in [13]. To obtain the statisti- 
cal accuracy of about 2% the sample of at least 150 electrons per crystal is 
needed. SND acquires such a sample daily when VEPP-2M operates in the 
center-of-mass energy 2E ~ 1 GeV. The average difference in calibration 
coefficients obtained using cosmic and e + e~ calibration procedures is about 
4%, while the e + e~ calibration improves the energy by 10%. For example, the 
energy resolution for 500 MeV photons improves from 5.5% to 5%. 



3.3 Energy and angular resolution 

The calorimeter energy resolution is determined mainly by the fluctuations of 
the energy losses in the passive material before (0.17Xo) and inside (0.17Xo) 
the calorimeter and leakage of shower energy through the calorimeter. The 
most probable value of the energy deposition for photons in the calorimeter is 
about 93% of their energy (Fig. 6). 
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In order to compensate for the shower energy losses in passive material and 
improve energy resolution the photon energy is calculated as: 

E = ai ■ Ei + a 2 ■ E 2 + a 3 ■ E 3 , (1) 



where E\ is energy deposition in the first and second layers of the central 
tower of the shower, E 2 is energy deposition in the first two layers outside 
the central tower, E 3 is energy deposition in the third layer, on are energy 
dependent coefficients. Here the tower is the three counters of the 1, 2 and 3 
layers with the same 9 and coordinates and the central tower corresponds 
to the shower center of gravity. 

The a, coefficients were determined from simulation of photons with energies 
from 50 to 700 MeV. For each photon energy the objective function: 

M = J2(E* - E k )\ (2) 
k 



was minimized over cej. Here E* is a known photon energy, Ek is the energy 
calculated using expression (1), k is a photon number. The energy depen- 
dences of ctj were approximated by the smooth curves. The approximation 
was done separately for the showers starting in different layers. As a result the 
calorimeter resolution was improved by 10%. 

The apparatus effects, such as nonuniformity of light collection efficiency over 
the crystal volume and electronics instability also affect energy resolution: 

<j e /E(%) = <ti(£) © a 2 (E) © <j 3 (E), (3) 



where cr 1 (E) is the energy resolution obtained using Monte Carlo simulation 
without effects mentioned above, cr 2 (E) is electronics instability and cali- 
bration accuracy contribution, <t 3 (E) is the contribution of nonuniformity of 
light collection over the crystal volume. For example, for 500 MeV photons 
a E /E = 5%, &i(E) = 3%, a 2 {E) = 1.2% and a 3 (E) = 3.8%. The depen- 
dence of the calorimeter energy resolution on photon energy (Fig. 7) can be 
approximated as: 

°b/E(%) = 4 /' 2% ■ (4) 



The distribution function of energy deposition in the SND calorimeter outside 
the cone with the angle 9 around the shower direction was obtained using 
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Monte Carlo simulation: 

E(9) = a-exp(-^fe/p), 



(5) 



It turned out that a and (3 parameters are practically independent of photon 
energy in a wide energy range 50 -j- 700 MeV. The method of estimation of 
photon angles based on this dependence was introduced in [14]. The depen- 
dence of the angular resolution on the photon energy shown in Fig. 8 can be 
approximated as: 

82° 

a = . ' © 0.63°. (6) 
y^GeV) 



The two-photon invariant mass distributions in — > 777 and — > 7i + 7i~7i° 
events (Figs. 9, 10), show clear peaks at tt° and 77 mesons masses. Invariant 
mass resolution is equal to 11 MeV for tt° and 25 MeV for 77. The kinematic 
fitting [15] improving angular and energy resolution of the detector. For exam- 
ple for — ■> 777, 7r°7 decays kinematic fit improves two-photon invariant mass 
resolution by a factor of 1.5 (Fig. 11). 



3.4 Particle identification 



The discrimination between electromagnetic and hadronic showers in the calorime- 
ter is usually based on difference in total energy deposition (Fig. 12). Multi- 
layer structure of the SND calorimeter provides additional means of particle 
identification based on differences in longitudinal energy deposition profiles. 
The distributions of energy deposition over layers for e ± and ^ is shown in 
Fig. 13. Two areas of pions concentration in this scatter plot correspond to 
nuclear interaction of pions and pure ionization losses in the first two layers of 
Nal(Tl). Utilizing differences in energy depositions for electrons and pions the 
special discrimination parameter was constructed. In the p(770) energy region 
it provides 99% selection efficiency for e + e~ — > 7r + 7r~, keeping contamination 
by e + e~ — > e + e~ events ~ 1% 

The distribution of energy depositions in the calorimeter layers is also used for 
Kl identification, for example in — > K$Kl decays, making feasible studies 
of rare K s decays. 

The e/n and 7/i^L separation parameters based on the differences in energy 
deposition profiles in transverse direction were also constructed. Their detailed 
description is given in [16]. 
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4 Tracking system 



4-1 Tracking system layout 



The tracking system (Fig. 14) consists of two cylindrical drift chambers and 
a cylindrical plastic scintillation counter (CSC) [17] between them. 

The counter length, inner diameter and thickness are 37, 12.6 and 0.5 cm 
respectively. In azimuth direction it is divided into 5 segments. The counter 
has a wavelength shifter fiber readout. It provides the time synchronization 
with beams collisions and produces signals for the first-level trigger. The time 
resolution of the counter is about 1.4 ns. 

The inner ("long") drift chamber (LDC) has the length of 40 cm, outer and 
inner diameters of 4 and 12 cm respectively. The corresponding dimensions of 
the outer ("short") drift chamber (SDC) are 25, 14 and 24 cm. 

Both chambers consist of 20 jet-type drift cells. Each cell has 5 gold-plated 
tungsten sense wires with a diameter of 20 /im. A ±300 fim staggering of 
the sense wires in azimuth direction resolves the left/right ambiguity. Field 
wires are made of bronze-coated titanium with a diameter of 100 /im. The 
total numbers of sense and field wires in each chamber are 100 and 260 re- 
spectively. The drift chambers flanges are made of fiberglass with a thickness 
of 10 mm. The wires are fixed using copper pins, precisely positioned in the 
flanges. The cylindrical fiberglass walls of both chambers have three layers 
of copper electrodes: field strips, parallel to the wires, and sense (cathode) 
strips, transverse to the wires, with electrodes for signal output. The effective 
thickness of tracking system is about 0.1X . Both chambers operate with a 
90%Ar + 10%CO 2 gas mixture. 

The circuitry of the electronics channel of the tracking system is shown in 
Fig. 15. The field wires and strips operate at different potentials in the range 
of 1-3 kV to provide the uniform drift field. The sense wires operating at 
ground potential are DC coupled to preamplifiers, while for cathode strips 
capacitive AC coupling is used. The amplified analog signals from each sense 
wire are transmitted to the amplitude- and time-to-digital converters (ADC 
and TDC). So, two amplitudes from both ends and a common timing are 
measured for each wire. 
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4-2 Calibration of the tracking system 



The calibration of the tracking system consists of several procedures to obtain 
the coefficients for conversion of measured times and amplitudes into the track 
coordinates and specific ionization losses (dE/dx). 

The sense wires are calibrated electronically using the generator. The gener- 
ator signals with varying amplitudes are switched between the preamplifiers 
on both ends of the wires. For each generator amplitude an average ADC 
count is calculated. These data provide calibration constants for ^-coordinate 
measurement by means of charge division method: 

, B L Q L - B R Q R 

: = l °- Ql + Qr ■ (7) 



Ql(r) — ^l(r)Ul(r) is the charge collected at the left (right) end of the wire, 
Ul(r) is the measured amplitude in ADC counts, A L ^ is the coefficient for 
conversion of the amplitude in ADC counts into units of electric charge, -Bl(r) 
is a conversion coefficient from charge into a normalized coordinate, l is the 
wire length. 

The ^-coordinate is also measured by means of cathode strip readout. To 
obtain a conversion coefficients from the ADC counts into the input signal 
amplitudes, the electronics channels are calibrated using generator. The track 
coordinate can be calculated then as: 

z = Z a + AZ, (8) 



where Z a is the coordinate of the strip with the maximum amplitude (central 
strip), AZ is the correction obtained from the ratio of the amplitudes at all 
fired strips to the amplitude at the central strip. AZ = f(q a -i, Q a , Qa+i), where 
a is the central strip number. Function / was obtained using experimental 
cosmic muon sample, under assumption of their uniform distribution in z. 

The small differences in wire lengths ^-coordinates of their centers contribute 
to the accuracy of z-coordinate measurements by charge division. To eliminate 
this contribution a calibration procedure based on analysis of e + e~ — > e + e~ 
events was performed. Cathode strips hit by one track were fitted and correc- 
tion coefficients for each wire obtained. The corrected z-coordinate reads: 

z = z w -(l + AL/L) - Az , (9) 



where z w is the track coordinate measured by charge division, AL/L is the 
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corresponding wire length correction, Azq is the longitudinal bias of the wire 
with respect to the center of the tracking system. 

Due to absorption in the gas mixture the average signal amplitude on a wire 
decreases by 1.5 times per 1 cm of the drift distance. The dependence of the 
average amplitude on drift distance can be approximated as: 

A(R) = A (l-R/l), (10) 

where R is a measured drift distance. Attenuation length I is obtained using 
e + e~ — > e + e~ events and A is the calibration constant. 

For the conversion of the drift time into the distance R(T) from track to the 
wire the r— calibration based on uniformity of distribution of e + e~ — > e + e~ 
events was performed. 

4-3 Drift chambers efficiency and resolution 

The efficiency of wires and cathode strips at low counting rates is close to 
100%. The drift chambers total counting rate is about 150 kHz for the beam 
currents I e + = I e - ~ 50 mA and luminosity of ~ 2 ■ 10 30 cm~ 2 s _1 . The pile up 
hits degrade accuracy of coordinates measurements. In order to increase the 
accuracy the hits distant from the track, are removed from the fit. Thus the 
final wire efficiency in track reconstruction is about 95%. The accuracies of 
coordinate measurements by drift times and charge division methods is about 
180 /im and 3.6 mm respectively. 

The tracking system angular resolution was measured using distributions in 
A6 and A0 acollinearity angles between tracks in the e + e~ — > e + e~ and 
e + e~ — > 7r + 7r~ events as a e = a Ae /\^2 and = where a is by 

definition the full width of the distribution at half-maximum divided by 2.36. 

The polar angle resolution o~q depends on the accuracy of z-coordinate mea- 
surement. The use of cathode strips together with sense wires provides 1.5 
times improvement of the resolution relatively to charge division only. Finally, 
oq = 1.7° for electrons and o~g = 1.9° for pions. The azimuth angle resolution 
o$ = 0.51° for electrons (Fig. 16) and = 0.54° for pions. The impact pa- 
rameter in r — plane resolution <tr = 0.5 mm. Small differences in angular 
resolutions for electrons and pions can be attributed to the different average 
ionization losses and 9 distributions for these particles. 

The achieved dE/dx resolution is about 30% and it is not sufficient for e/n 
separation, but 7r ± jK ± separation in the energy region of 0-meson production 
is possible (Fig. 17). 
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5 Muon system 



The muon system provides suppression of cosmic background and identifi- 
cation of the muons produced in the beams collisions at energies above 450 
MeV. 

5.1 Muon system layout 

The muon system (Figs. 1, 2) consists of 16 (200 x 40 x 1 cm 3 ) scintillation 
counters and two layers of streamer tubes. Each counter consists of two sheets 
of plastics scintillator glued together. Light is collected from two sides of the 
scintillator via the acrylic light guides to the pair of PM tubes with photo- 
cathode diameter of 40mm. Each counter is wrapped in aluminized mylar and 
put in a thin steel container. 

The barrel part of the streamer tubes system consists of 14 modules with 40 cm 
width, two end-cap modules have the width of 80 cm. The module length is 
2 m. Each module consists of 16 tubes arranged in two layers. The tubes are 
made of 300 fxm stainless steel and have a diameter of 4 cm. Along the tube 
axis a 100 /im gold-plated molybdenum wire is strung. Streamer tubes are 
filled with a 75%Ar + 25%n-pentane gas mixture, which provides operation in 
a streamer mode at wire potential ~ 4000 V. 

The electronics channel of the muon system [18] is shown in Fig. 18. The analog 
signals from PMTs are carried to the discriminators where the amplitude and 
timing signals are formed, and then transmitted to ADC and TDC. The signals 
from opposite ends of the streamer tubes are applied to the time expander, 
where the signal with a duration proportional to their time difference is formed 
for further digitization by TDC. The signals from the tubes and counters are 
also used in the first-level trigger. 

5.2 Muon system calibration 

Z-coordinate of the particle track is obtained using the time difference between 
the signals from the opposite ends of streamer tubes : 



where / is the wire effective length, known with the accuracy of about 2%, AT 
is the time difference between the signals from wire ends expressed in TDC 




C — AT 

C 



(11) 



10 



counts, C is the calibration coefficient, obtained using AT distribution for 
cosmic muons. Stability of the electronics was checked daily using calibration 
generator. 

To obtain the time of scintillation counter hit the following expression is used: 

t — K-T + At, (12) 

where T is a TDC count, K is a coefficient for TDC count to time conver- 
sion, obtained from cosmic muons spectrum, At is a correction accounting 
for differences in the individual cables lengths and PM tubes delay times. 
At is obtained as the time difference between the hits in muon scintillation 
counter and tracking system scintillation counter, using e + e~ — > events 
and cosmic muons. Cosmic muons are also used for the energy calibration of 
the scintillation counters (average energy losses are 2 MeV/cm). The energy 
deposition in the counter E = k ■ U, where U is the signal amplitude in ADC 
counts, k is the coefficient of conversion from ADC counts to the units of 
MeV. To consider the amplitude dependence on the Z-coordinate of the muon 
track due to the light absorption in scintillator the muon track is fitted using 
streamer tubes. 



5.3 Efficiency and time resolution of the muon system 

The muon system veto signal provide the 50 times cosmic background sup- 
pression in actual experimental conditions. 

A longitudinal coordinate of the track is obtained from the tube number. If 
only one tube is fired, then resolution a = dj \f\2 = 1.2 cm, where d = 4 cm is 
the tube diameter. When the tubes in both layers are hit the resolution is a = 
0.6 cm determined by the shift between layers. The Z-coordinate resolution 
o~z = 3.3 cm. 

To obtain the time of scintillation counters hit with respect to beams crossing, 
special corrections for light propagation time were applied. The corrections 
take into account the effective velocity of light in the scintillator (~ 15 ns/cm) 
and dependence of the counter firing time on the signal amplitude due to fixed 
discriminator threshold. After corrections the time resolution of the counters 
is 1 ns (Fig. 19). 



6 Data acquisition system 
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6. 1 Electronics 



The SND data acquisition system (DAQ) [19] (Fig. 20), based on KLUKVA 
electronics standard [20] developed at BINP. 

Analog signals from the detector elements come to the front-end electronics lo- 
cated on SND systems. Then the signals are carried via screened twisted pairs 
to digitizing electronics (ADCs and TDCs). The signals from the calorime- 
ter and cathode strips are shaped by shaping amplifiers before digitization. 
The digitizers and shapers reside in KLUKVA crates. One crate holds 16 data 
conversion modules, one first-level trigger (FLT) interface module (IFLT) and 
readout processor module (RP). The data conversion modules provide data, 
which are transferred via KLUKVA crate bus to IFLT modules for further use 
in FLT. 

Two types of digitizers are used: TDCs with "common stop" and ADCs with 
"common start". The first type [21] provides the drift time measurement. 
Signals from sense wires start 250 MHz counters and FLT signal stops drift 
time digitizing. If FLT signal does not arrive within 1 /is the counters reset. 
The second type [22] is used for digitization of signals from the calorimeter 
and DC cathode strips. The digitization starts by FLT signal and lasts 80 /is. 
The digitizers (about 3000 channels) occupy 13 KLUKVA crates. 

After the digitization the data are extracted from the ADC and TDC modules 
by the RP in 100 ns. An RP memory contains the pedestals of all channels 
placed in the crate. RP compares the data with pedestals and performs zeros 
suppression. As a result average event size decreases down to about lkB. 
The RP internal memory keeps digitized data for two events, providing with 
the internal registers of ADCs and TDCs 3-level derandomizing of data flow. 
The time needed for event digitization and processing is about 200 /is. The 
data from RP are transmitted to the ON-LINE computer via special CAMAC 
modules [20]. 

6.2 First-level trigger 

The fast analog signals of energy deposition in calorimeter towers and logical 
signals from wire hits in drift chambers are collected by IFLT modules. The 
IFLTs and analog and logical summation modules form the signals for spe- 
cialized trigger modules: Calorimeter Logic, Track Logic, Layers Logic, and 
Energy Thresholds, producing 48 different FLT components (Fig. 21). The 
input information for Calorimeter Logic [23] is 20 signals from sectors and 
8 signals from rings, produced by the logical summation of the hit towers 
with the energy deposition higher than 25 MeV. The sectors signals go to 



12 



the address bus of 1MB static RAM, where the look-up table is stored. The 
information from rings from another 32 kB RAM is added. So, according to 
the RAM contents, the components of the FLT depending on the number of 
hit towers and their relative positions are produced. FLT components can be 
easily changed by RAM reprogramming. 

The Energy Thresholds module contains 10 discriminators with programmable 
threshold. As an input information the analog signals of total energy deposi- 
tion and energy deposition at "large" angles are used. 

A Track Logic [24] is constructed using the programmable logic devices. As the 
input information the following logical signals are used: 100 signal lines from 
LDC wires, 20 signals from the first two layers of SDC and 20 signals from the 
calorimeter sectors. Hit wires patterns in LDC corresponding to tracks from 
interaction point are looked for. Then the track continuation in SDC and 
calorimeter is checked. Resulting signals contain information on the number 
of found tracks and their relative positions. The Layer Logic is based on static 
RAM and as input data uses 10 logical OR signals from LDC and SDC layers. 
It produces three components, each being the various combination of input 
signals, according to the RAM contents. 

The FLT logic is implemented as a pipeline working at 16 MHz clock rate, 
which is the collider beam circulation frequency. Such FLT organization pro- 
vides zero "dead" time, while the FLT decision latency is about 800 ns. 

The main source of background are electromagnetic showers produced in the 
collider magnets and lenses by stray particles from the beam. The FLT back- 
ground suppression is based on the following differences in the properties of 
physical and background events: 

(1) the total energy deposition in the calorimeter for physical events is larger 
than for background events; 

(2) hits in physical events form two or more compact clusters in the calorime- 
ter; 

(3) particles in physical events originate from the collider interaction region; 

(4) in the background events the main part of the total energy deposition in 
the calorimeter is localized at "small" angles. 

The most often used FLT components are: 

(1) "OR" of all calorimeter towers used as a master signal for FLT Mask 
Modules with a time resolution better than 5 ns, while the time between 
collisions is 60 ns; 

(2) two well separated towers hit; 

(3) two collinear towers hit; 

(4) hit tower at an angle larger than 54° with respect to the beam; 
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(5) two hit towers with azimuth angle between them larger than 90 degrees; 

(6) the total energy deposition and the energy deposition at "large" angle is 
larger than certain threshold; 

(7) LDC track, LDC track continued in SDC, and LDC track continued in 
SDC and calorimeter; 

(8) two separate tracks in LDC; 

(9) the muon system veto. 

Final FLT decision is produced by 10 Mask Modules, each taking 48 FLT com- 
ponents as inputs. It is possible to define 10 different trigger types by writing 
appropriate masks into the Mask Modules. Four trigger types select events 
containing charged particles ("charged" triggers), other four select events with 
neutral particles only ("neutral" triggers), one trigger selects cosmic muons 
for muon system calibration and the last selects the e + e~ — > e + e~ events for 
on-line luminosity measurement. For all "neutral" triggers the muon system 
veto is enabled. 

In average experimental conditions the calorimeter and drift chambers count- 
ing rates are about 80 kHz and 150 kHz respectively. The FLT reduces the 
trigger rate to a level of ~ 80 Hz. 

6.3 Data Taking 

Data from RP RAMs are read into the VAXserver 3300 computer via the 
CAMAC interface modules at a rate of 2 ms/event. The computer RAM is 
used for data flow derandomization. Then the data are transmitted to the 
VAXstation-4000/60 workstation, where events are packed in a special SND 
format [26]. For packed events the partial reconstruction (see section 7) is 
performed. The reconstructed events are processed by the software 3-rd level 
trigger which suppresses beam background and cosmic events, decreasing the 
event rate by a factor of 2. In addition it counts e + e~ — > e + e~ and 77 events 
for on-line luminosity monitoring. The total time required for partial recon- 
struction and 3-rd level trigger decision is ~ 10 ms/event. Average DAQ dead 
time is about 10%. 

Events selected by the software trigger are written on file on the local disk. 
After each experimental run files are automatically saved to 4 GB tapes by 
the data management system ART [19,27]. 

During the run several concurrent processes collect information about detector 
subsystems, DAQ, and collider performance. Real-time access to these data is 
provided by information management system IMAN [28]. The special process 
monitors this information and rises alarm if hardware or software problems 
are detected. 
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Stability of SND subsystems is monitored using precision generator signals to 
test all electronics channels. This procedure takes about 20 minutes and is 
performed daily. 



7 Data processing 

Large volume of data collected during experimental runs requires quick and 
adaptable processing. The data processing is based on CERN software such 
as HBO OK [30], HIGZ [31], PAW [32], MINUIT [33], as well as programs 
developed in BINP: 

(1) GIST [29] provides common environment for data processing, user inter- 
face with other codes, data storage in format compatible with HBO OK 
package for further analysis using PAW; 

(2) COCHA [26] provides the data input/output and storage. The data rep- 
resentation is based on entity - relationship model; 

(3) event reconstruction code; 

(4) code for events visualization based on HIGZ graphics package; 

(5) ART [27] provides transparent access to data stored on tapes and disks. 

Processing of experimental data consists of several steps: 

(1) SND subsystems are calibrated using e + e~ — > e + e~ and fi + fi~ events; 

(2) events are reconstructed and parameters of particles are calculated; ac- 
cording to these parameters the events divided into ~ 20 different physical 
classes; 

(3) for each experimental run the e + e~ — > e + e~ and 77 events are counted 
for the integrated luminosity determination, efficiencies of SND systems 
are obtained; collected information about SND and collider performance 
is used in the analysis, for example, to determine the beam energy more 
precisely; 

(4) the selected events are used for the physical analysis. 

Event reconstruction is performed in a following sequence. First step is a 
search for separated clusters in the calorimeter. Then the track reconstruction 
in drift chambers is performed. Tracks are linked to the calorimeter clusters. 
The clusters with energy depositions of more than 20 MeV not linked to 
tracks in DC are considered as photons. For charged tracks no requirements 
on energy depositions in the calorimeter are imposed. Finally parameters of 
reconstructed particles are calculated: angles, energies, dE/dx losses. 

The computer system (Fig. 22) used for reconstruction and analysis of exper- 
imental data includes several desktop Pentium PCs running Linux 2.0 and a 
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Sun Enterprise 450 server under Solaris 2.6 operating system. All computers 
are connected via 100 Mb Ethernet to a local network. 

The server provides common software and user disk space for desktop PCs, 
along with batch system and network services. It also runs event reconstruction 
code at an average speed of about 150 events per second. 



8 Monte Carlo simulation 

The Monte Carlo simulation of the SND is based on UNIMOD [35] package. 
SND geometrical model comprises about 10000 distinct volumes. The primary 
generated particles are tracked through the detector media taking into account 
the following effects: ionization losses, multiple scattering, bremsstrahlung of 
electrons and positrons, Compton effect, e + e~ pair production by photons, 
photo-effect, unstable particles decays, interaction of stopped particles, nu- 
clear interaction of hadrons. After that the signals produced in each detector 
element are simulated. To provide the adaptable account of experimental con- 
ditions: the electronics noise, signals pile up, the actual time and amplitude 
resolution of electronics channels. Broken channels are also taken into account 
during reconstruction of Monte Carlo events. 

Comparison of experimental and simulated distributions are shown in Figs. 23, 24. 



9 Conclusions 

The design and performance of the SND detector operating at VEPP-2M 
e + e~ collider since 1995 is described. During this time the total integrated 
luminosity of ~ 25 pb~ l in the center-of-mass energy range 2Eq = 0.4 -h 1.4 
GeV was accumulated. New rare radiative decays <j) — > 7r°7r°7 [36], r]7r 0r y [37] 
and OZI and G-parity suppressed decay — > cj7t° [38] were observed. The 
(f) — > r} >r ) [39] decay existence was confirmed. Many other decays and hadronic 
production cross sections were measured with higher accuracy. More integrated 
luminosity is to be accumulated in the region of p and u> mesons providing 
data for studies of their rare decays. 
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Fig. 1. SND detector, section along the beams: (1) beam pipe, (2) drift chambers, (3) 
scintillation counter, (4) light guides, (5) PMTs, (6) Nal(Tl) crystals, (7) vacuum 
phototriodes, (8) iron absorber, (9) streamer tubes, (10) 1 cm iron plates, (11) 
scintillation counters, (12) and (13) elements of collider magnetic system. 




Fig. 2. SND detector, section across the beams: (1) beam pipe, (2) drift chambers, 
(3) scintillation counter, (4) Nal(Tl) crystals, (5) vacuum phototriodes, (6) iron 
absorber, (7) streamer tubes, (8) scintillation system. 
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Fig. 3. Nal(Tl) crystals layout inside the calorimeter: (1) Nal(Tl) crystals, (2) vac- 
uum phototriodes, (3) aluminum supporting hemispheres. 



FLT 




Fig. 4. Electronics channel of the SND calorimeter: (Nal) Nal(Tl) scintillator, (VPT) 
vacuum phototriode, (CSA) charge-sensitive preamplifier, (ADC) analog to digi- 
tal converter, (Gen) calibration generator, (SHA) shaping amplifier, (ATT) com- 
puter-controlled attenuator, (FLT) first-level trigger 
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Fig. 5. Energy spectra for photons with energy Eq = 500 MeV after cosmic calibra- 
tion. E is a measured energy. 




Fig. 6. Dependence of the most probable energy deposition on the photon energy. 
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Fig. 7. Dependence of the calorimeter energy resolution on the photon energy, 
oe/E(%) = 4.2%/ \/E( GeV). E - photon energy, ge/E - energy resolution ob- 
tained using e + e~ — > 77 (dots) and e + e~ — > e + e~7 (circles) reactions. 
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Fig. 8. Dependence of the angular resolution on the photon energy, 
0-0 = 0.82%/£( GeV) 0.63°. E - photon energy 



22 




Fig. 9. Two-photon invariant mass distribution in experimental 4> — > 777 events. 
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Fig. 10. Two-photon invariant mass distribution in experimental 4> — > 7r + 7r~7r' 
events. 
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Fig. 11. Spectra of the invariant masses of photon pairs in e + e — > 777 events 
before (dots) and after (line) kinematic fit. 
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Fig. 12. Energy deposition spectra for 385 MeV and ir^ 
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Fig. 13. Energy deposition in calorimeter first two layers vs total energy deposition 
for 385 MeV e ± and vr ± . E - total energy deposition, E\-\- E% - energy deposition 
in the first two layers 




Fig. 14. SND coordinate system view across the beams: (1) beam pipe, (2) scintil- 
lation counter, (3) field shaping wires, (4) anode wires, (5) field shaping strips 
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Fig. 15. Electronics channel of Drift chambers (DC): (1) DC gas volume, (2) end 
cap, (3) pins for wire attachment, (4) field shaping wires, (5) sensitive wires, (6) 
field shaping strips, (7) sensitive strips, (HV) potential from high voltage divider, 
(CSA) charge sensitive amplifier, (PA) preamplifier, (A24) 24 channel amplitude to 
digital converter, (T2A) time and two amplitude to digital converters 




Fig. 16. The azimuthal acollinearity angle distribution for e + e — > e + e events 
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Fig. 17. dE/dx distribution for events with charged particles in the center-of-mass 
energy range 2Eq ~ 1 GeV 
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Fig. 18. Muon system electronics: (AC) scintillation counters, (Tubes) streamer 
tube, (D) discriminator, (EXPAN) time expander, (Trigger) first-level trigger, (X) 
logical summator, (TDC) time to digital converter, (ADC) analog to digital con- 
verter 
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Fig. 19. Time distribution of scintillation counters signals with respect to the beams 
collision time in e + e~ — > n + n~ events at 2Eq ~ 1 GeV. The cosmic muons back- 
ground is also seen as a uniform contribution. 
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Fig. 20. SND data acquisition system: KLUKVA crates, (FLT) first-level trigger, 
(IFLT) FLT interface modules, (Logic) modules of FLT logic and decision, (CA- 
MAC) CAM AC crates 
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Fig. 21. First-level trigger (FLT) structure scheme: (IFLT) FLT interface modules, 
(FADC) ADC "total energy deposition", (A24) ADC, (S) - logical summator 
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Fig. 22. The computing system. 
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Fig. 24. Angle V> between pion pairs from the Ks — > vr + 7r decay 
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